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SUMMARY 

Equations are developed to show that, for one-dimensional two-group 
neutron-distribution calculations, the properties of a neutron -ref lector 
region can be combined in three parameters that effectively represent the 
action of the reflector. The use of these parameters and boundary con- 
ditions based on them for single and multizone reflectors is considered 
and from the results of computations, the following applications of the 
parameters are shown: 

(1) Use of the parameters in connection with reactor core proper- 
ties for comparison of reflector configurations 

(2) Direct use of parameters as circuit constants in electrical 
analog 

( 3 ) Use of boundary conditions based on reflector parameters to 
obtain numerical solutions of two-group equations in a neutron -producing 
region. 


INTRODUCTION 

The close analogy between the physical concepts and equations for 
the diffusion of neutrons in a reactor and the flow of current in an 
electrical transmission line suggests the application of already 
developed electrical engineering techniques for handling problems con- 
nected with reactor design. The use of lumped electrical circuit 
elements to approximate a continuous transmission line has been success- 
fully employed to simulate the diffusion of neutrons by electrical analog 
simulators (reference l). In the analogy, the neutron fluxes and the 
diffusion currents are simulated, respectively, by the electrical 
potentials and the currents of the transmission line. The following 
analysis and examples carried out at the NACA Lewis laboratory describe an 
additional transmission-line concept that can be applied to the one- 
dimensional solution of the neutron -group equations. 

The neutron diffusion, absorption, and moderating properties of the 
material in a passive reactor zone of finite extent give rise to a p^ti- 
cular relation between neutron fluxes and diffusion current at the interface 
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between that medium and an adjacent reactor zone. The diffusion of 
neutrons into the passive zone can be treated in identical manner to the 
flow of ciirrent into an analogous electrical transmission line connected 
to an unspecified sending circuit. The transmission line is ordinarily 
characterized by a definite lengthy terminating impedance, characteristic 
impedance, and attenuation constant.. In order to specify its reaction on 
the sending circuit, the line is replaced by an equivalent sending -end 
impedance that behaves like the line with respect to its reaction on any 
other electric circuit. The neutron reflector, when described by a similar 
set of. properties, can be treated in the same analytical manner as the 
transmission line. In the two-group approximation the only processes 
occurring in a passive reactor zone are the reflection of fast and thermal 
neutrons and the moderation and reflection of fast neutrons as thermal 
neutrons. A minimum of three reflector parameters will, therefore, describe 
the action of the passive zone. The use of transmission-line-like proper- 
ties gives this set of three parameters and introduces the following advan- 
tages and simplifications. 

(1) Reflectors or general passive zones with structures based on 
usual engineering requirements can be directly compared by a consider- 
ation of these three parameters for each in connection with the require- 
ments imposed by the reactor core. (A simplified example presented for 
water and beryllium reflectors demonstrates the comparison of reflector 
parameters. ) 

(2) The definition of a neutron-diffusion impedance for a reflector 
region is immediately applicable to the nuclear -reactor simulator. Here, 
the reflector network is actually replaced by a single resistor or combined 
resistor and current source. The simplification leads to more rapid 
operation of the simulator and permits the use of network sections usually 
assigned to the reflector to give a more detailed measure of the neutron 
distribution in the active region of the reactor. 

(3) The effective replacement of a reflector by suitable properties 
applied at the core-reflector interface suggests the use of these proper- • 
ties in numerical calculations to define new boundary conditions applied 
only to the neutron flux in the reactor core (appendix A) . ^ This simplified 
boundary-condition concept for a simply reflected reactor leads to a set 

of two boundary equations in place of the usual four equations. The 2X2 
criticality determinant derived from the boundary equations is, however, 
of approximately the same complexity as the ordinary criticality determi- 
nant reduced in order to a 2X2 determinant. The simplified boundary 
conditions are especially useful when applied to a multiregion reflector 
(appendix B) . There is then no increase in the complexity of the criti- 
cality determinant. 
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ANALYSIS 

Neutron-diffusion properties ♦ - The definition of characteristic 
diffusion properties for a neutron-reflector material is based on the 
neutron -balance equations applied to the reflector region of a nuclear - 
reactor segment with its center at the origin of the veiriable r . Con- 
sider first the equation for the fast-neutron group: 



- ckt“ Za,f ^f 


2 

"3 ® 


= 0 


( 1 ) 


The constant oj determines the size of the reactor segment considered; 
m is an integer, 1, 2, or 0, respectively, depending on whether the 
reactor geometry is cylindrical, spherical, or rectangular. The 

factor cur^ then serves to make equation (l) a neutron-conservation 
equation applicable to the entire reflector segment and permits a com- 
parision with the analogous electrical-transmission-line equation with 
a direct nitmerlcal equivalence of terms . The nuclear constants are 
defined in the list of symbols (appendix C), and the buckling constant B 
applies only to reactor configurations having finite linear dimensions 
in other than the r -direction. Equation (l) is similar in form to the 
equation for the potential on the pictured nonuniform electric trans- 
mission line with distributed resistance and leakage. 



The line of length (r^-r^.) has a series resistance per unit length 
Rq and a shunt conductance per unit length A poten- 
tial is applied at the sending end of the line at r^. The 

terminating impedance on the line is zero or infinite, respectively, 
for the requirements of =0 or of S0j.(r|j)/Sr = 0. These 

requirements correspond to those for a finite neutron reflector bounded 
by (l) a nonreflecting region, or (2) a perfect reflector or identical 
image reflector region and source of neutrons. Equation (l) transforms 
to the transmission-line equation by the definition of the following 
quantities : 



(2a) 



(2b) 



4 


KSCA RM E52H01 


The transformed equation (l) is then 


+ m 

Ro^f dr2 r dr 


“•rf r oi - n 

T — - ^ ^0,f ?*f - 0 


(3a) 


or its equivalent 




dr 


+ “ __ 

2 r dr 


" ®0,f = 0 


(3b) 


Tyo combinations of the quantities defined in equations (2) appear 
in equation (3b) eind its solution. Therefore let 


ry f-i U • J 

^o,f = fH- 
^o,f 


(4a) 


Of^ = Ro,f ^o,f 


(4b) 


Because of their transmission-line analogs, these quantities are defined 
as a fast-neutron characteristic impedance ^O.f and a fast -neutron 
attenuation constant Of . 

The neutron-balance equation for the thermal group corresponding to 
equation (l) is 


03T' 


.m 


^tr,th ^^^th m "^^th 

C _L 


dr 


F" r ~d?“ “ ^a,th 0th - ojr“ 


^r,th 


0th 


= - oxr“ Pth^af 0f 


(5) 


This equation also transfonns into a transmission-line equation that is 
equivalent term-by-term with the thermal -group equation. 


r“ f ^^J^th ^ m ^^th 


^0,th V dr^ 


+ “ 


r dr 


• ^ ^0,th 0th = - r™ g^h^f 0f 


(6a) 


or 


— ^ • °0,th ^0,th 0th = - ^,th Sth,f 0f 

cir^ i ui 


r dr 


(6b) 
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with the substitution of the following quantities: 


^0,th 



®0,th = “2:a,th 


. ‘“\r,th 


(7a) 

(7b) 


Sth,f = "^th £a,f 


(7c) 


The thermal -neutron source indicated on the right side of equation (6a) 
is analogous to a current introduced per unit length in an electrical 
transmission line. It is therefore helpful to define further nuclear 
parameters appearing in equations (6) and in their solution that correspond 
to transmission-line parameters. Let 


Z ^ = 
0,th 


Rp,th 

®0,th 


(8a) 


a 


th 


2 _ 


= R, 


0,th 


(8b) 


These thermal parameters correspond to the fast-neutron parameters defined 
in equations (4) . 

Reflector solutions . - Equations (3b) and (6b) with m = 0, 1, 
or 2 are to be solved satisfying either of the following two 
boundary conditions depending on the desired reactor configuration: 

(1) The flux vanishes at the outer boTindary of the reflec- 
tor (0(r^j) = 0) 

(2) The gradient of the flux is zero at the origin (0'(O) = O) 

The general solution of equation (3b) satisfying the appropriate 
boundary conditions is 

0jf(r) = Cf <t(cuj.r) (9) 

where is a constant determined from the known flxox at the reflector - 

core interface r^, and ®(ofr) is a generalized function the specific 

value of which is given in table I. 

The net neutron diffusion current Ip(r) in the r-direction at 
any point in the reflector is 
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If(r) = - ^ = - r” ^ f(c^r) 


’’o.f 


0,t 


where 


The ratio 


0f(r) 


$ '(or) s -- — - ♦(otr) 
d(or) 


is determined entirely by reflector properties and 


can be defined as the fast-neutron impedance of the reflector (r) : 

( 10 ) 


Z,(r) 

^ 4>’((Xfr) 


Specific values of Z|.(r) appear in table I. 
Then 

0f(^c) 


If(rc) 


Zf(rc) 


( 11 ) 


defines a new boundary condition at r = r^ for the neutron flux in the 
adjacent region. 

The solution of equation (6b), 

^^$*^th(^^) m ‘^0th(’^) 


dr 


+ 

2 r 


dr “ “"th ^th^’^^ ^ ■ ^0,th ®th,f 


for the thermal neutron flux is 




°f^ ■ °4;h^ 


( 12 ) 


where 4>((xthr) is the generalized homogeneous solution of (6b) given in 
table I for the reflectors considered. 
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The net thermal diffusion current Ith(^^) 'tlie r-direction is 

* p dr 

^0,th ^ 


or 


Ith(r) = - 4”(othr) + ^th^f^f ^ $,(^^ 3 ,) 


4 -v, ' ' ' ?. 2 

Equation (13) can he put into a more convenient form with the 
definition of a thermal -neutron impedance Zth(r) : 


(13) 


^th(^) 


^0,th ^(°th^) 
r"> ^’(athr) 


(14) 


which . corresponds to Z^(r) defined in equation (lO) . 

With the use of this thermal -neutron impedance, equation (l3) 
becomes 


Ith(r) 


^th(r) 

Zth(’^) 


^0,th gth,f ^f 

2 2 

Of - octh^ 


4>(ctfr) 

Zth(r) 



or equiveilently at the reflector boundary r„. 


Ith(rc) - Pth,f(rc) If (re) 

th' c’ 


or 


Ith(’^c) = 


_ ^th^^c) ^th,f (^c^ 


0f(^c) 


where 


^th,f(’^c) = 


R, 


0,f 


2-f (rc) 


R 


0,th ^th(’^c) 


"gth,f ^0,th ) 

- o^h^ 


(15a) 


(15b) 


(16) 


Equations (ll) and (l5) state that the net fast-neutron and thermal- 
neutron diffusion currents into the reflector are linear functions of the 
fast- and thermal -neutron fluxes at the reflector boundary. 
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Equations (ll) and (l5b) can be used as a new set of bo\mdary 
conditions for the two -group equations of the core. 

Two-group criticality equations . - The fast- and thermal -neutron 
diffusion ciirrents are given by equations (ll) and (l5b) for a particular 
reactor segment. Division of equations (ll) and (l5b) by axT^n results 
in two equations giving the neutron diffusion currents per unit area. 
These equations applied at the core -reflector interface as boundary 
conditions to the general solution of the two-group core diffusion equa- 
tions (appendix A) 

0^(r) = (P^r) + D 02 (pgr) (A7) 


0tl,(r) = A 0^ (p^r) + B Sg ©2 (Pg^) 
lead to the criticality determinant 


L 3 Sr £cr“Zf (r)J 

^r,th „ S®i(3ir) 

■ Sr 

. [ ' 


[- 


^tr,f ®2^P2^^^ 


Sr asr™Zf(r)_ 


03r“Zj(r)y 


0l(Pir) 


^r,th „ ^2(^2^) 




Pth(r) 


oa'“Zf(r) 


- 02(^2^) 


(A8) 


= 0 (17) 


r=Tp 


The solution of this determinant can be obtained by varying any of 
several variables: (l) the reflector thickness in Z(r), (2) the radius 

of the core, or (3) the fissionable-material concentration in S 
the p . 


The ratio of the parameters in equations (A7) and (A8) is 


\r,f 

S0^(p^r) 

0l(Pir) ' 

3 

6r 

(crm Z^(r) 

r^tr,f 

302 (Pgr) 

^ 02 (P2^) ■ 

L 3 

hr 

oor™ Zf(r)_ 


(18) 


This ratio leads to knowledge of the spatial variation of 0f(r) and 
0th(r) in the core . 
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APPLICATION OP REFLECTOR PARAMETERS 

. Comparison of reflector properties . - A group of reflectors 
satisfactory from an engineering and nuclear standpoint can be further 
compared by a knowledge of Pth,f > Zf , and Zth^ which give a descrip- 
tion of the over -all behavior of the reflector. These quantities are 
particularly useful for a comparison of reflector properties in the 
case of a multireflector reactor where the problem is difficvilt to 
solve by ordinary two -group methods. 

The following properties of ZfCr^-r^), Z^hCr^-rc), and Pth^fC^b-rc) 
facilitate a choice of reflector for a particular core structure. 

(1) A low Zf(r^-rg) will depress the fast flux in the core leading 

to a large neutron current into the reflector region, whereas a high 
Zf(rb-rc) tends to confine the neutrons to the core where they will 

slow down . 

(2) A large Z^j^(r^-rg) tends to prevent the diffusion, and hence 
the loss, of thermal neutrons from the core. 

(3) The quantity Pth,f ® Pleasure of the moderating power 
of a reflector, a large value of Pth,f(^b"^c^ indicating considerable 
slowing down, resulting in more thermal neutrons returning to the core. 

A demonstration of these properties for simple water and beryllium 
reflectors of various thicknesses is given in figures 1 and 2j spherical 
geometry and a core radius of 30 centimeters are assimed. A water 
reflector thickness of 10 centimeters gives values of Zf and Z^j^ 

95.1 and 99 percent of the respective saturation values while the cor- 
responding values for 20 centimeters of beryllium are 98.1 percent and 

89.2 percent. The saturation value of ^ for water, 0.286, means 

that not more than 28.6 percent of the fast neutrons can slow down to 
return to the core. The graphs therefore indicate that little is gained 
by choosing a water reflector thickness of greater than 10 centimeters 
or a beryllium thickness of greater than 20 centimeters for this 
particular reactor configuration. Although the beryllium reflector has 
a lower value of Z^^ the water reflector, the beryllium reflector 

would be more desirable on the basis of its higher value of Zf and Pfh f 
providing reflector thickness is no problem. 

Simulator operation using reflector parameters . - A typical reactor 
criticality determination was performed on the NACA nuclear reactor 
simulator for a simply reflected cylindrical reactor in order to establish 
a procedure for the direct use of reflector parameters. For the fast- 
neutron group, the reflector was replaced by a fast-neutron resistance 
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equal to Zf(rg) that terminated the electric netvork for the active 

zone. Measurements made on the network then included that of the cxirrent 
throu^ this terminating resistance. The thermal -neutron network for 
the reactor core was then terminated hy a resistance equal to Z^jj(rg), 

the thermal -neutron diffusion impedance. In addition, a current Pth,f 

multiplied hy that measured through the reflector fast-neutron resistance 
was introduced at the terminal point of the core network. 

The over-all reflector parameters give the following advantages in 
simulator operation: 

(1) The network sections ordinarily assigned to the reflector can 
he used to give more detailed information of the neutron distribution 
in the reactor core. 

(2) The reduction in time required to set one resistor value and 
one current value in contrast to a large number of adjustments required 
for a reflector network greatly simplified the operation of the simulator. 
This feature is especially useful for reactor calculations involving the 
adjustment of core parameters hy iterative processes. 

(3) Various reflectors can he employed in calculations hy the simple 
readjustment of one resistor value for each neutron group. 

(4) If a detailed variation in flux for a multizone reflector is 
required, the analytical expression with boundary values measured on the 
simulator gives more conplete information than can he found with most 
practical reflector networks. 

Numerical calculation example . - In order to determine computational 
time, a numerical calculation based on the criticality determinant 
(equation (l7)) was made for a spherical reactor having a single reflector 
and compared with the computational time involved in solving the two -group 
core-reflector equations . 

Equation (17) has in it several possible variables: (l) the radius 

of the core r^, in the parameters Z^, P-th,f> reflector 

structure in the parameters Zf, Z^j^, and (3) the fissionable 

material concentration in the parameters S, p, and ®. 

Any of these quantities, for example, the fissionable material 
concentration, may be used as a variable j however, rc was made to vary 
as with the usual computational method in order to obtain a direct 
comparison. Because of the dependence of Z and P^h^f ^c^ 

2x2 determinant is approximately equal in complexity to the reduced 
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fonn of the 4X4 determinant and it was found that there is no loss or 
gain in computational ease or time between the two iterative methods . 

For a single core having internal and external passive zones, the 
concept of reflector impedances and a f (r) leads to a 

4X4 determinant) and there is no increase in the order of the 
determinant due to the addition of more passive zones . The conven- 
tional solution to the same problem requires an 8X 8 determinant with 
4 additional rows and col\amns for each additional passive zone. Thus, 
a considerable saving in computational effort and time should result 
from the use of the reflector impedance method over the conventional 
method of solution of this particular problem. 


CONCLUSIONS 

A neutron-reflector region can be effectively represented by only 
three parameters for two-group calculations in one dimension. For the 
fast-neutron group a single quantity, the fast-neutron impedance, is 
necessary, and for the thermal group, a thermal -neutron impedance and 
current parameter are required. The simplified boundary conditions 
dependent on these parameters have the following advantages: 

1. Consideration of the three reflector parameters in connection 
with a proposed core structure permits a direct comparison of the 
effectiveness of reflectors that is especially useful for miultizone 
reflectors . 

2. Reflector parameters are applicable as actual circuit elements 
on an' electrical network simulator for the group diffusion model. 

3. Boundary conditions based on the reflector parameters can be 
applied to flux calculations- in neutron -producing regions without further 
consideration of the neutron flux in the reflector regions . 


Lewis Fli^t Propulsion Laboratory 

National Advisory Comimittee for Aeronautics 
Cleveland, Ohio 
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APPENDIX A 


DERIVATION OF CRITIC ALITY EQUATIONS FOR REACTOR-PLUS -REFLECTOR 
CONFIGURATION BASED ON TWO-GEOUP TBEORY AND 
TRANSMISSION-LINE THEORY 


The resiilts of the .work on the finite transmission line can he 
applied to the two-group equations for a reactor -plus -reflector 
configuration. Essentially, the reflector equations and the four core- 
ref lector interface boundary conditions are replaced by two boundary 
conditions j this, in effect, reduces the ordinary 4x 4 determinant 
solution of the diffusion equations to a 2X2 determinant solution. 


The one -dimensional two-group equations for a reactor core are 
(reference 2 ) 

v/ |!lf (r) - i«j(r) - 0^(r) 


+ k. . Z 0.. (r) =0 

th a,th '^th' ’ 


(Al) 


and 



^tr,th t>2 


® - ^a,th * I’th = 0 


(A2) 

p " 

where B represents the buckling constant for directions normal to rj 
it applies only to the solutions involving rectangular and cylindrical 
geometry and does not apply to a spherical geometry. 


The two boiindary conditions replacing the two reflector equations 
and the four core -reflector interface boundary conditions are 


hr,f ^0f(r) 0f(^^) 

■ 3 dr - 

and 

_ ^tr,th *^^th(^) ^ 0th ^th,f(^) v 
3 dr “ oir’Ji Z^j^(r) ' osrm Zf. 


(A3) 


(A4) 
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Equations (Al) and (A2) must satisfy the following relation; 


0^(r) = 5^f(r) 


(AS) 




(A6) 


Substitution of equations (A6) and (AS) in equations (Al) and (A2) 
will give two equations in the buckling constant (reference 2). 

This constamt p^ may have either of two values: one value positive 



and one value negative 



The general solution of equations (AS) and (A6) will be given by 


0f(r) = A ©1 (Pir) + D ©gOgi*) (A7) 

0^j^(r) = A Sj0^(p^r) + D (pgr) (A8) 

where ©^^(Pj^r) and © 2 (P 2 r) represent the value of the solution for 
^f(r) and 0^j^(r) corresponding to the two values P]^^ and Pg^^ 
respectively. The constants S^^ and S 2 represent the ratios of 


©th/Qf^ the components of 
respectively . The values 

the fluxes for the values of and 

of © 2 ^ and ©2 are 

^2^ 

s — - 

^tr,f 


^a,f 

[1 - kf(l-Ptij| 

(A9) 

^1 - 

2 .. 
th "a,th 

^a,th 


^r,f 

(?2^-b2) + 

"a,f ___ 

[1 - kf(l-Pth)] 

(AlO) 

^2 - - 

^^th ^a,th 

^th^a,th 


The values for 0]^(P]^r) and © 2 (P 2 r) are listed in the following 
table for m = 0, 1, 2: 
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Vr^ 

0^(p^r)> p^^>0 

esCPjr); 9/<0 

m = 0 

dr^ 

cosh (p^r) 

COS (P2I*) 

m = 1 

d2 Id 
^2 r dr 

Iq (3i^) 

Jq (P2^) 

m = 2 

d^ 2 ^ 

^2 r dr 

sinh (P]^r)/r 

sin (p2r)/r 


Substitution of the boundary conditions ((A3), (A4)) into equa- 
tions (A7) and (A8) at r = will give 


^tr,f 

3 ' 


S 0i(Pir) 5 

A -• ■ + D ^ 


hr 


hv 


■] 


air™Zj(r) 


[a 0^(p3^r) + D ©gCPg^)] 

(All) 

o^riaZth(r) 

+ D Sg ©2(^2^)] - ^ ®2(P2^)] (A12) 

Simplification of equation (All) and (A12) results in two equations in 
the parameters A and D: 


^tr,th 

3 


[as^ 


^ ®l(Pl^^) ^ ^ ^ ®2(P 2^) 

Sr 2 




[• 


>^tr,f ^ ®i(3ir) _ 0i(Pir) 

3 Sr asr“-Z^(r) 


]*•[- 


V,f ^ ®2(P2^)\ ®2(P2^) 




?ii!2lLlD = o 

(0r®E^(r) 


(A13) 
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^tr,th „ ^ Qi(Pir) / Si Pth,f(^) \ ^ 

3 1 Sr ■ W“Zth(r) 'asr^Cr)] 


r) 


g S 02(p2^^) _/ Sg 
^ \aar”'z^j^(r) asr®Z^(r) 


ttr™Z^(r)/ 


02(32^) 


D = 0 


(A14) 

There will be a nontrivial solution to equations (A13) and (A14) only if 
the determinant of the coefficients vanishes. Thus, the general 
criticality equation is 


[ 

■[ 


^tr,f ^ 

3 hr CDr^^(r) 




S2 


hr 


®2(P2") 


^tr,_f ^ _j_ ®2(P2^) 




(r) 


>^tr,th s. ^ 


bv 


yaxr®Z.(.j^(r) ccir°^£.(r)/ 2 

- f-ii - Wll)0i(Pir)|= 0 

Va3r®Z|.jj(r) u3r“Zf.(r)/ J 


(A15) 


Depending upon the type of problem under consideration, the reflector 
thickness, the radius of the core r^, or the fissionable material con- 
centration may be varied to achieve criticality. 

The ratio of the two constants in equations (A7) and (A8) is 


f\tr,f ^ ®i( 3 ir) © 3 ^( 3 j^r) 


B 1 

_ 3 Br cjor^^(r)J 

A “ 

^■tr,f ^ 


.3 Sr mr“Zf(r)_ 


If the ratio d/a is known, the spatial variation of flux in the 
core may be plotted as a function of the distance r. The value of the 
flux in the core will then be known to an arbitrary constant. With the 
values of the fast- and thermal -neutron flxoxes known at the core-reflector 
interface, the reflector fluxes may be calculated by using equations (9) 
and (12) with the constants adjusted to give the proper flux values at 
the core-reflector interface. 
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APPENDIX B 


TWO-GROUP PARAMETERS FOR MULTIPLE PASSIVE ZONES 

The reflector consisting of 'two regions as shown schematically in 
the following figure can he represented hy three parameters similar to 
those for a simple one-region reflector; 


(Region a) 

(Region b) 

^c 

i 


Region h, extending from r^ to the boundary is immediately 

replaced by a set of parameters as described in the body of the report 
or, if it is a multizone region, by a set previously determined by the 
following process. These parameters establish necessary boundary 
conditions for the region-a neutron flux solutions (r^, ^ r ^ r^^) euad 

lead to reflector parameters that describe the effect of the entire 
reflector combination# 

Fast-group parameters . - The general one-dimensional solution for 
the fast-neutron flux in region a has the form 

= Af[ilri(afr) + Df 

where Aj and D^ are constants to be determined by the applied boundary 
conditions, and the functions i|r 2 ^(a^r) and \lf 2 (ct^r) are listed in the 
accompanying table for three coordinate systems: 


Coordinate 

system 

ti(ctr) 

^ 1^2 (dr) 

Rectangular 

slnh (ctr) 

cosh (oo?) 

Cylindrical 

Kq (or) 

Iq 

Spherical 

sinh (ctr) 

cosh (ctr) 

r 

r 


At r^, the boundary of region a, the neutron current is 


I(r ) = 

f i Ro,f dr 


r “ A 
^1 ^f 

Zo,f 


• (o^ri) + Df \lf2 • (o^ri)] (B2) 
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This current is related to the flux at the boimdary r. by the 
in^jedance ^ ^ 

_ _ Zo,f 

If(ri) " 

Boundary condition (B3) permits evalxiation of the constant 

jj ^ _ Zp,f ti’(gfr^) 

^ ^0,f 'l^2(“'f’^i) + ^i“ ^f,h i2’(°f^i) 

The fast-group flvix and current solutions obtained with permit 

the definition of a fast-neutron ii^edance for the reflector combination 
according to 


tlCofTi) + Df TlTgCo^r^) 


tl ' (o^i) + Df tg ’ (ofri) 


(B3) 


Z 


f,c 


^fC^c) 

If(^c) 


The value of this impedance is 


^f,c = 


^0,f 

, m 


ti(oLprc) + Df ^|f2(gfr^,) 

♦i'(Vc> * “f ♦a'Vcl 


(B5) 


Therma3.-group parameters . - The thermal-group equation for region a 
is 


J^th - “th^ 5^th = - Ro,th Sth,f Af [^^iCafr) + Df i|f2(ofr)| (B6) 

The general solution of equation (B6) is separated into two parts con- 
sis ting of a complementary function 

5^th,h = Cfi, [il^iKhr) + &th 1'2Kh’^)] (B7) 

and a part made up of the. particular integral and an arbitrary multiple 
of the complementary function 
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where ®th constants . This separation of the 

solution leads to equations that simplify the definition of parameters 
for the double reflector. 

The bo\indary conditions for the two parts of the thermal -group 
solution are 


^th,h(^c) ^ i^th(^c) 


T fy> \ ^5^th,h 

Hh,h^^i^ = " R 

* ^0,th ^ 


^ 0th,h(^j) 


"th,b 


(B9a) 

(B9b) 


0th,p (^c) = 0 


(BlOa) 


^ R0,th ^ 


- Pth,f,b If(ri) (BlOb) 

^th,b " ' 


Evaluation of bo\mdary condition (B9b) gives 


■'th,b ” 


^0,th 

^i 


m 


^l(“th^i) + ^th ’•^E^^th^i) 


from vdiich the constant 6. 


th 


'I'l ’ (Wi) + ^th '^ 2 ' (“th’^i) 

is evalxiated 


^th = ■ 


^0,th ^l(°th^j) •*• ^i° ^th,b ^l'(°th^j) 


L^O,th ’*^2(°'th^i^ + ^i”" ^th,b *2'(^h’^i^J 


(Bll) 


The thermal -neutron diffusion in5>edance for the combined reflector is 
defined iii terms of the homogeneous part of the flux solution in 
region a. The current of thermal neutrons into the reflector 




is written 




^ \ th^ c . 

th,c 


thjp^’^c^ 


(B12) 
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The homogeneous part of equation (B12) is then solved for ^ with 

condition (B9a) 




^0,th * \h 

'•' 2 ' Kh''c> 


(B13) 


The imdetermined constants in the particular flux solution (B8) 
are fixed by the boundary conditions (BIO) . Equation (BlOa) gives 


^h = 


^O.th Sth.f ^f(’^c) 


(of^-«th^) [ti(o^h’^c) + ®th ^2(°4:h^c)] 


(B14) 


This value of A^j^ permits the evaluation of boimdary condition (BlOb), 
giving 

V, ^ ^ ’•"I'Kh^) + Dth '1^2’Kh^i) ^ , 

-h + n 0c (^c) 

l^lKh^c) + ®th 'l'2(‘=4ih’^c) 


= I 0f(rc) 


tl(°4;h^l) + Dth ^^2(«th^j) 
^l('"th’"c^ + ^h ^2(‘"th^c^ 


(B15) 


^i°^ ^0,th 6th ,>f 

^0,th 


(B16a) 


I ^ ^0,th 6th, f 

^th,b^“'f^'“'th^^ 


(B16b) 


^0,th ^th,f 

^th, ^(of2_a^j^2) 


^th,f ^0,f ^ ^th,f ,b ^i(°-f^i) + Pf tgCa^r^) 

^f^bC'^^-Oth^) ti(oi-rc) + Bf teCofi-c) 


(B16c) 
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The constant then is obtained directly from equation (B15) : 

_ ^ ^ ’I'lC'Hhrc) + ^ ’t'l'Kh^i) 


I>th = 


“ ^2(°'th’'c^ - ^ ^ ^2’Kh^^ 

This permits a complete evaluation of p(rg): 


^th,p(^c^ “ 


(B17) 


It.. ^ 

h ^c 

^l'(Wc) + ^th ^2'Kh^c)' 

\ 

^0,f Sth,f 


_^l(“th^c) + ®th ^2(°th’^c) - 



r M^c) 

(B18) 


In the form eiaployed for the single reflector region this is equivalent 
to 


Ith,p(^c) = - 0f (’^c) (B19) 

^f,c 

The three constants ^f,c> B-th f ^ appearing in equa- 

tions (B5) , (b 13) , and (B19) , respectively, then serve to represent 
all properties of the double reflector for two-group calculations. 
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APPENDIX C 
SYMBOLS 


The following symbols are used in this report; 
A,C,D,6 undetermined constants 


8th,f 

I(r) 

k 

m 


buckling constant for directions normal to r (applies to 
solutions involving rectangular or cylindrical geometry) 

neutron conductance defined by equation (2b) 

neutron transconductance defined by equation (Vc) 

neutron diffusion current at point r 

neutron multiplication constant 

0 , 1 , 2 


^th,,f 

Pth 

RO 

r 


reflection parameter defined by equation (16) 
neutron resonance escape probability 
neutron diffusion resistance defined by equation (2a) 
distance variable 

position of reflector outer boundary 


position of core-reflector interface 


S 


Z 


Z 


0 


a 


3 


2 


position of boimdary between two reflector regions 
ratio of flux components ©th/®f 
reflector impedance 

characteristic impedance of reflector material 

attenuation constant of reflector material 
buckling constant for core (r -direction) 
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Ko(a^b)Ao(°^b) 

general function 
neutron transport mean free path 

macroscopic neutron absorption cross section (includes slowing 
in fast group) 

general function 

neutron fliax at point r 

general function 

reactor segment size factor 

m d. 
dr2^^^ 

Subscripts : 

b parameter for region b 

c parameter for combined reflector 

f fast neutron group 

h indicates solution of homogeneous equation 

\ 

p indicates particular solution 

th thermal neutron group 
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Figure 1. - Fast and thermal impedances for beryllium and water reflectors. 
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